home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / cgecon.z / cgecon
Text File  |  1996-03-14  |  3KB  |  133 lines
  1.  
  2.  
  3.  
  4. CCCCGGGGEEEECCCCOOOONNNN((((3333FFFF))))                                                          CCCCGGGGEEEECCCCOOOONNNN((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      CGECON - estimate the reciprocal of the condition number of a general
  10.      complex matrix A, in either the 1-norm or the infinity-norm, using the LU
  11.      factorization computed by CGETRF
  12.  
  13. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  14.      SUBROUTINE CGECON( NORM, N, A, LDA, ANORM, RCOND, WORK, RWORK, INFO )
  15.  
  16.          CHARACTER      NORM
  17.  
  18.          INTEGER        INFO, LDA, N
  19.  
  20.          REAL           ANORM, RCOND
  21.  
  22.          REAL           RWORK( * )
  23.  
  24.          COMPLEX        A( LDA, * ), WORK( * )
  25.  
  26. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  27.      CGECON estimates the reciprocal of the condition number of a general
  28.      complex matrix A, in either the 1-norm or the infinity-norm, using the LU
  29.      factorization computed by CGETRF.
  30.  
  31.      An estimate is obtained for norm(inv(A)), and the reciprocal of the
  32.      condition number is computed as
  33.         RCOND = 1 / ( norm(A) * norm(inv(A)) ).
  34.  
  35.  
  36. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  37.      NORM    (input) CHARACTER*1
  38.              Specifies whether the 1-norm condition number or the infinity-
  39.              norm condition number is required:
  40.              = '1' or 'O':  1-norm;
  41.              = 'I':         Infinity-norm.
  42.  
  43.      N       (input) INTEGER
  44.              The order of the matrix A.  N >= 0.
  45.  
  46.      A       (input) COMPLEX array, dimension (LDA,N)
  47.              The factors L and U from the factorization A = P*L*U as computed
  48.              by CGETRF.
  49.  
  50.      LDA     (input) INTEGER
  51.              The leading dimension of the array A.  LDA >= max(1,N).
  52.  
  53.      ANORM   (input) REAL
  54.              If NORM = '1' or 'O', the 1-norm of the original matrix A.  If
  55.              NORM = 'I', the infinity-norm of the original matrix A.
  56.  
  57.      RCOND   (output) REAL
  58.              The reciprocal of the condition number of the matrix A, computed
  59.              as RCOND = 1/(norm(A) * norm(inv(A))).
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. CCCCGGGGEEEECCCCOOOONNNN((((3333FFFF))))                                                          CCCCGGGGEEEECCCCOOOONNNN((((3333FFFF))))
  71.  
  72.  
  73.  
  74.      WORK    (workspace) COMPLEX array, dimension (2*N)
  75.  
  76.      RWORK   (workspace) REAL array, dimension (2*N)
  77.  
  78.      INFO    (output) INTEGER
  79.              = 0:  successful exit
  80.              < 0:  if INFO = -i, the i-th argument had an illegal value
  81.  
  82.  
  83.  
  84.  
  85.  
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.